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We show that the two-channel Anderson lattice model leads to the de- 
velopment of an unconventional superconductivity out of a metallic non 
Fermi-liquid phase. It is characterized by a composite order parameter 
comprising of a local spin or orbital degree of freedom bound to triplet 
Cooper pairs with an isotropic and a nearest neighbour form factor. The 
gap function is non analytic and odd in frequency, and a pseudo-gap devel- 
ops in the conduction electron density of states which vanishes as \uj\ close 

to U! = 0. 

PACS numbers: 75.20.Hr, 75.30.Mb, 71.27.+a 



1. Introduction 

Heavy Fermion [1] (HF) superconductivity (SC) has drawn much atten- 
tion since the discovery of superconductivity in CeCu2Si2 [2] which is likely 
characterized by an anisotropic order parameter with symmetry yet to be 
determined. It became apparent over the last decade that almost all HF 
materials are unstable with respect to magnetic or superconducting phase 
transitions, which either compete with each other or can even coexist as 
found in uranium based materials. Moreover, superconductivity in UBei3 
develops out of an incoherent metallic phase with a high resistance at T c of 
~ WO/iQcm. The two-channel Anderson lattice model, 

+ E aV {c\ aa X% + h.c} , (1) 

is believed to serve as a possible descriptions of the electronic properties of 
UBei3 or Pr 3+ compounds with cubic symmetry, eg is the band dispersion of 
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the two conduction bands a coupling to the localized doublets via hybridiza- 

(i) 

tion V. X a g are the usual Hubbard operators describing local /-states at 
each lattice side i. One doublet is assoziated with a quadrupolar charge 
distribution (X aiOC ) and one carried a spin (X a>a ). Its paramagnetic phase 
is characterized by a large residual resistivity and entropy, and ill defined 
electronic quasi-particles. Fermi liquid physics is restored by cooperative 
ordering or applied magnetic field [3]. 



2. Composite order parameter 

It was observed that the local two-particle-particle irreducible vertex 
has al/w singularity and is odd with respect to the incoming and outgoing 
frequencies [4]. This is related to the fact that correlations are induced 
through consecutive hybridization processes on the /-shell due to the strong 
on-site Coulomb repulsion. Therefore, odd-frequency superconductivity [5] 
is favoured in this model. Assuming an order parameter with T± symmetry, 
the order parameter must have spin singlet, channel single (SsCs) symmetry 
or spin triplet, channel triplet (StCt) symmetry. Only for the latter 
stable solution was found. The order parameter in the StCt sector 

Oij = sign(£ CT - E a ) [(siPj'*) - foif a >] , (2) 

correlates a local /-spin (quadrupolar moment) (rj) - i, j = x,y, z - com- 
ponent with a itinerate Cooper pair P s,t (P t,s ) in the SsCt (StCs) sector: 

ps,t = ^Y, s ^^)%%m-k) (3) 

k 

pt ' S = jf E S(k)^(k)iz y %i j(-k) ■ (4) 
k 

The form factor S(k) transforms as Ti, and tp(k) is a bi-spinor in spin and 
channel space. Oij is invariant with respect to exchanging the orbital and 
the magnetic doublet \a) and \a) and a simultaneous particle-hole trans- 
formation of the conduction electrons. Hence, the occurrence of supercon- 
ductivity is independent of the character of the local ground state of the 
local ion. We investigate the superconducting phase within the dynamical 
mean field theory (DMFT) by self-consistently calculating the anomalous 
self-energy under the presents of a Cooper-pair field. 



3. DMFT in the superconduction phase 

The Nambu Green function is an 8 x 8 matrix in spin and orbital 
space. Assuming no directional coupling between spin and orbital degrees 
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of freedom, the anomalous 4x4 self-energy matrix may be written as 
g(z, k)a 2 z 2 n s dn c T_ where n s and n c are constant unity vectors in spin 

and channel space, and the amplitude function g(z, k) has to be odd in fre- 
quency. This reduces the full problem to a standard 2x2 size for g(z, k) 
and the diagonal self-energy S. In order to derive DMFT equations with 
a purely local self-energy matrix S (z), the anomalous self-energy is re- 
stricted to isotropic pairs, e. g. g(z,k) = g{z). G_ c (z) denotes the medium 
matrix in which the effective impurity is embedded. It is related to a gen- 
eralized Anderson width matrix through A(z) = V 2 g_ 2 G_ c (z) g_ 2 and has 
normal and anomalous components describing quasi-particle propagation 
and pair-creation and annihilation, respectively. The band electron self- 

energy = T_ I + GT is determined by the local T-matrix T_ whose 

diagonal elements are given by the local quasi-particle scattering matrices 
T_ n (z) = V 2 Gf(z), and its anomalous contribution stems from Cooper pairs 
scattered on the local shells. 

As an impurity solver for effective site problem of the DMFT, the NCA 
[6] was extended to the superconducting phase. The additional terms AE 7 (z) 
to the self-consistency equation for the local ionic propagators are generated 
by the anomalous media in the SC phase. 

/oo roc 
/ d ede'p f (e)p f (e)f(-e') (5) 
-OO J —OD 

■P a {z + e)P a (z + e - e')P a {z - e') 

/ dede'p f (e)p f (e)f(-e') (6) 

-oo J — OO 

■P a (z - e)P a (z + e' - e)P a (z + e') 

where p/(e) = 9m[G]i2(e — i5)/ir is the spectral function of the anoma- 
lous media. The anomalous local T-matrix is also obtained by the same 
additional diagram of the generating functional and reads 

/oo p e —(3z 

^dep f (e)f(-e) f c ^-^P a (z) (7) 

■P a (z + iuj n )P a (z + iu n - e)P a {z - e) , 

where Zf is the partition function of the effective site, and the contour C 
encircles all singularites of the integral kernel. These coupled double integral 
equations are numerically solved using similar techniques as developed for 
the Post-NCA equations [7]. 

We calculated the self-consistent DMFT solution of the spectral func- 
tions close to T c which are displayed in Fig. 1. The linear frequency de- 
pendency of the quasi-particle and anomalous Green function very close to 
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Fig. 1. Spectral functions of the quasi-particle and the anomalous Green function 
in DMFT(NCA) in the superconducting phase for 1 — T/T c <C 1 in the vicinity of 
u> — fj, = 0. The dotted curve shows a fit with a* \lu\. Parameters: p(e) = exp(— [(e — 
fi)/W] 2 )/{y^W) conduction band DOS, band width W = 10A, E a -E a = -2A, 
A = V 2 Trp(v), 

lo = 0, similar to d-wave superconductivity, is clearly visible. The opening 
of a pseudo-gap indicates a decrease of the energy in the superconducting 
phase, and, hence, its thermodynamical stability. Details of the filling de- 
pendence and the phase diagram are beyond the scope of this paper and 
will be published elsewhere [4]. 
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